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Abstract. We show new data from the 64 Ni+ 124 Sn and 58 Ni+ 112 Sn reactions studied in 
direct kinematics with the CHIMERA detector at INFN-LNS and compared with the reverse 
kinematics reactions at the same incident beam energy (35 A MeV ). Analyzing the data with 
the method of relative velocity correlations, fragments coming from statistical decay of an 
excited projectile-like (PLF) or target-like (TLF) fragments are discriminated from the ones 
coming from dynamical emission in the early stages of the reaction. By comparing data of 
the reverse kinematics experiment with a stochastic mean field (SMF) + GEMINI calculations 
our results show that observables from neck fragmentation mechanism add valuable constraints 
on the density dependence of symmetry energy. An indication is found for a moderately stiff 
symmetry energy potential term of EOS. 



1. Introduction 

Heavy ion collisions are a powerful tool to study the nuclear properties at different conditions of 
density, temperature and isospin asymmetry. For example, in the Fermi energy regime (10-100 



A MeV ), it is possible to explore from very low barionic densities (p/po < 0.2) clustered matter 
PQ to densities in proximity to the saturation value (po = 0.16/m~ 3 ) [2j. Conversely, in the 
relativistic energy regime (E/A > 100 A MeV ) it is potentially possible to access densities 
up to 2 — 2.5 times the saturation value in a short timescale Hj. For these reasons it is 
not surprising that heavy ion collisions, in particular with projectiles and targets with large 
isospin asymmetries, have been widely used to probe the density dependence of the symmetry 
term of the nuclear Equation of State (EOS), that is a key ingredient for the dynamical model 
calculations of heavy ion collisions and astrophysical predictions [5J. Different experimental 
observables have been used to constraint the density dependence of the symmetry energy in 
heavy-ion collisions: isospin diffusion and equilibration [6], neutron to proton ratio [7|, light 
charged particles transverse collective flow [5] , ratio of fragments yields and isoscaling [9] , heavy 
residues production in semi-central collisions |10j . isospin migration in the low density "neck" 
region [llj . These different observables from heavy ion reactions at low densities are generally 
consistents within uncertainties |12j . but results can be strongly dependent from comparison 
with microscopic calculations simulating the nuclear dynamics. Indeed a weak overlap exists so 
far between constraints from heavy ion collisions and the estimate of the symmetry energy from 
astrophysical observations on neutron stars [131 [T2] . For these reasons it is important to add 
new observables, increasing accuracies of the existing ones and to improve the model simulations 
in order to reduce the relative large actual uncertainties. 

In this contribution we present experimental data for the 64 Ni+ 124 Sn and 58 Ni+ 112 Sn 
reactions studied in direct kinematics with the CHIMERA detector at the same beam incident 
energy (35 A MeV ) of the previously studied experiment in reverse kinematics [141 [TT] . The 
ensemble of data of the two experiments collects a unique set of information on the midrapidity 
"neck" fragmentation mechanism in semi-peripheral dissipative collisions. We show that the 
Intermediate Mass Fragments (IMF, 3<Z<20) midrapidity emission presents many experimental 
properties (like the N/Z isospin asymmetry enhancement) that in transport models calculations 
are attributed to the formation, in the early stage of the reaction, of a low density region 
("neck") connecting projectile-like and target-like fragments. These properties can be linked to 
reaction dynamics if they are correlated with the timescale evolution of the nuclear reactions. 
In this contest we have compared experimental data for the reverse kinematics experiment with 
a stochastic mean field calculation (SMF) in order to get a parametrization for the potential 
symmetry energy term of EOS. 

2. Experimental details 

The direct kinematic experiment has been performed at the Catania INFN-LNS superconducting 
cyclotron where beams of 64 Ni and 58 Ni at 35 A. MeV of 0.5 enA intensity impinged on self- 
supporting 124 Sn and 112 Sn thin targets whose thickness are respectively 187 fig/ cm 2 and 
268 fig /cm 2 . Reaction products were detected with the 4ir detector CHIMERA in its full 
configuration, constituted by 1192 Si-CsI(Tl) telescopes arranged in 35 rings around the beam 
axis. Data acquisition was triggered when at least three silicon detectors were fired. This 
experiment was the first to use on a large number of detectors the CHIMERA silicon pulse 
shape upgrade [15], giving the possibility to charge identify the particles that are stopped in 
silicon detectors. In the reverse kinematics experiment only the forward part of the detector 
(688 telescopes covering the angular range between 1° and 30°) was used. Details of this latter 
experiment setup are described elsewhere (Til [TT]. 

Fig. [JJ shows an example of the charge identification obtained for light fragments (Z<20) 
using the pulse shape discrimination (PSD) technique for a silicon detector. The identification 
threshold is around 4 A. MeV for light charges (Z=6). The use of PSD technique was very 
important in this direct kinematics experiment because a part of slow IMFs, produced in the 
mid-rapidity region, fall below the threshold for AE — E identification. We have selected almost 
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Figure 1. Rise-time of impinging particles vs. particles energy bi-dimensional plot for particle 
stopping in a silicon detector at 13.75° for the reaction 64 Ni+ 124 Sn. The inset shows the quality 
of charge identification obtained. 
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Figure 2. Atomic number of fragments Z as a function of their parallel velocity along the beam 
axis V par for the selected ternary events. PLF and TLF regions are indicated by labels. 



complete events where the total charge is 45 < Ztot <80 and the parallel momentum of the 
colliding system is at least 60% of the total one. Semipheripheral collisions were selected gating 
on the total charged particle multiplicity M<7. Following the same methods used in Refs. 
|144 [TTj in order to select ternary events in the final state, we considered a subset of events 
where the total charge of the three biggest fragment Z(l)+Z(2)+Z(3) > 45 and their momentum 
p(l)+p(2)+p(3) > O.Q Pbeam . 

Fig. [2] shows for the selected events the 2-dimensional plot of fragments parallel velocity 
along the beam axis as a function of their charge for the 64 Ni+ Sn reaction. It is interesting to 



compare this plot with the same obtained for the reverse kinematics experiment (Fig. 2 of Ref. 
|14j). Here the slow- moving and massive TLF fragments are produced around 1 cm/ns velocity 
and all particles belonging to this region (mainly stopping in the first stage silicon detectors) 
are identified in mass with Time-of-Flight techniques. The complete angular coverage at least 
up to 90° becomes essential to study events in which projectile-like and target-like fragments 
are simultaneously present in the same events for semi-peripheral collisions. We clearly see also 
the two regions respectively of PLF (Ni-like) with velocities around 7.5 cm/ns and of IMFs in 
the intermediate velocity region between the two main partners of the reaction. 

3. Results 

In order to disentangle dynamically and statistically emitted fragments and to look for timescale 
of fragment formation, the three biggest fragments of each event were sorted according to the 
decreasing value of their parallel velocity, following the method described in [11]. The three 
particles, labeled as PLF, IMF (3<Z<20) or TLF, depending upon their respective velocity were 
analyzed to check the correct attribution event by event; finally the fragment-fragment relative 
velocities Vrel(PLF, IMF) and Vrel(TLF, IMF) were calculated and are reported in Fig. [3j 
The relative velocities are normalized to the one corresponding to the Coulomb repulsion as 
given by the Viola systematics [16]. Fig. [3]shows the correlations between the two relative veloc- 
ities r 1 =V REL (PLF,IMF)/V V ioia(PLF,IMF) and r 2 =V REL (TLF, IMF)/V Viola (TLF, IMF) 
for the IMFs charges Z=4,6,10,14. As compared with a similar plot in the Reverse experiment 
(Ref. [2], fig. 7) we note now that we can populate with similar efficiency both the regions 
along the r\=\ axis (whose yield is dominated by sequential emission from PLF) and r 2 =l axis 
(whose yield is mainly due to sequential emission from TLF); values of r\ and r 2 simultaneously 
larger than unity indicate a prompt ternary division (dynamical origin). We can observe that 
heavier fragments (as Z=14 in the figure) are originated mainly by the break-up or fission of 





V REL /V VIOLA (PLF,IMF) 



Figure 3. For the Ni+ Sn reaction, correlations between relative velocities Vrel/V\ 
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of the three biggest fragments in the event for Z/mf=4,6,10,14. The points refers to a 
monodimensional Coulomb trajectory calculation (see text). 



the heavy partner, i.e. the target-like residue in our case and they lie along the T2=l axis, while 
light fragments are concentrated along the diagonal, indicating their prevailing dynamic origin. 
In Fig. [3] a timescale calibration was done, as in Ref. [2] using a three-body collinear Coulomb 
trajectory calculation. The inner points along the diagonal correspond to the shortest timescales 
(40-60 fm/c), corresponding to IMFs predominantly emitted from the dynamically expanding 
neck region formed at midrapidity, between the projectile-like and target-like primary fragments. 
This method permits to explore and disentangle different stages of the dynamical evolution of 
the system, that normally are mixed together when looking, for example in a source recognition 
analysis, at Galilean invariant Vji — V± contour diagrams. 

In order to study the alignment properties of midvelocity fragments, we have evaluated the 
angle 9prox- As shown in the inset of Fig. 4d), if the IMF had its origin from a PLF break-up, 
Oprox is the angle between the (PLF-IMF center of mass)-TLF relative velocity axis, and the 
PLF-IMF break-up axis (relative velocity between PLF and IMF oriented from the light to the 
heavy fragment). This definition mainly differs from the ones presented in Refs. |18t ll9| because 
it requires the explicit detection of a TLF and PLF fragments in the same event. cos(9ppox)=^ 
indicates a complete alignment with the IMF emitted in the backward hemisphere respect to the 
PLF and cos(9ppox) < indicates the emission in the forward hemisphere respect to the PLF. 
Of course the IMF emission could be attributed also to a TLF break-up. As a first approximation, 
we have calculated cos(9ppox) requiring the condition Vrel(PLF, IMF) /Vyioia <l-6. This 
condition removes most of IMFs emitted by TLF |17j . A more complete analysis in the direct 
kinematics experiment will give the possibility to extend and complete these results considering 
the contribution due to both the TLF and PLF break-up. 

Figdjshows the cos(9ppox) angular distribution for Z=4 and Z=8 IMFs for the two reactions 
under study in the direct kinematic reaction (a,b left panels) and in the reverse kinematic 
reaction (c,d right panels) respectively. The yields are normalized to the respective area. For 
a statistical emission the cos(9ppox) distribution is expected to show a forward-backward 
symmetry around cos (Oppox) = 0. We note that the distributions are peaked at cos(9ppox) ~ 1 
indicating a strong anisotropy favoring the backward emission respect to the forward one in a 
strict aligned configuration along the TLF-PLF separation axis. In Ref. [11] has been shown in 
details how, adding a kinematic condition on the fragments degree of alignment to a condition 
that selects the greatest deviations from Viola systematics in the relative velocities correlation 
plots (as illustrated in Fig. [3]), it is possible to careful disentangle the pattern of dynamically 
emitted fragments (characterized by short emission times and by a strong degree of alignment 
in cos (Oppox)), respect to statistically emitted fragments. For these two classes of events it 
is now possible to introduce a further observable, i.e. the fragments isotopic composition N/Z. 
This is the object of the last section. 

4. Isotopic distributions and comparisons with SMF model 

Fig. [5h) shows the < N/Z > as a function of the IMFs atomic number Z for the reaction 
124 Sn+ 64 Ni. Fragments statistically emitted in the PLF forward hemisphere have been selected 
by using the condition cos(9ppox) < 0. The relative points are shown as solid squares in 
Fig. Eh). Solid circles shows the < N/Z > for dynamically emitted IMFs. Events for 
these particles are selected by imposing that cos(9ppox) > 0.8 (highest enhancement for 
backward emission) and selecting events near the diagonal in the VpEL/Vvioia(PLF, IMF) — 
VpEL/Vvioia(TLF, IMF) relative velocities correlation plot. We note that the < N/Z > ratio 
is systematically larger for dynamically emitted particles respect to the statistically emitted 
ones. In Fig. [5b) we have reported (solid circles), for the same reaction, the correlation between 
cos(9ppox) an d < N/Z > for all fragments with charges between 5 < Z < 8. We observe 
an increase of the < N/Z > at values of cos(9ppox) approaching to 1, corresponding to the 
highest degree of alignment. A similar result has been found recently, in a different data analysis 
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Figure 4. a),b) panels: cos(9ppox) angular distribution calculated for a PLF (Ni-like) break- 
up and production of Z=4 and Z=8 IMFs for the reactions 64 Ni + 124 Sn (red histogram) and 
58 Ni+ ii2g n ^ blue higto 

gram); c), d) panels: same distributions for the for the reverse kinematic 
reactions 124 Sn+ 64 Ni (red histogram) and 112 Sn+ 58 Ni (blue histogram) for a PLF (Sn-like) 
break-up. The insert in Fig. 4d) gives a sketch of the Oprox definition. 



contest, for the 124 Xe+ 124 ' 112 Sn system [19] . 

The data, for the inverse kinematics neutron rich reaction 124 Sn+ 64 Ni, were compared with 
a transport theory using the stochastic mean field model (SMF) based on Boltzmann-Norheim- 
Vlasov (BNV) equation |20U21| . The potential part of the symmetry energy is taken into account 
using two different parametrizations as a function of density named asy-stiff and asy-soft. The 
first one linearly increases with the density while the second one exhibits a weak variation around 
the saturation density po- The slope parameter of the symmetry energy, defined as 



L = 3p 



dE sym (p) 



dp 



(1) 

P=P0 



is in the current calculation around 80 MeV for the asy-stiff and 25 MeV for the asy-soft choice. 
The statistical code GEMINI [22] is used as a second step de-excitation phase applied to the 
SMF primary hot fragments. 

In Figs. [5)3, c) the SMF + GEMINI calculations are plotted as hatched area histograms for 
dynamically emitted fragments. The hatched zone indicates the error bars in the calculation. 
As can be observed in Fig. [5fc>) the asy-stiff parametrization (blue hatched area) produces more 
neutron rich fragments respect to the asy-soft choice and the difference persists after the GEMINI 
secondary-decay stage for Z<7. The asy-stiff parametrization matches the experimental data 
fairly well. This is confirmed in Fig. [5b) where the asy-stiff parametrization better reproduces 
the < N/Z > enhancement observed for values of cos(9ppox) > 0.9. In the calculations the 
dynamical emitted fragments production yields become near to zero below cos(6ppox) <0.7, 



thus theoretical < N/Z > values can not be calculated for cos(9prox) below «0.7. In Fig. [5b) 
SMF calculations for primary fragments and asy-stiff parametrization are also shown (empty 
circles symbol). 

We have checked if the SMF + GEMINI calculation with asy-stiff parametrization reproduces 
other basic experimental observables. As an example, we show in Fig. [6] the experimental 
longitudinal velocity distribution respect to beam axis (in the cm. system) of Z=4 
dynamically emitted fragments (empty histogram). Velocity distributions (yellow histograms) 
of experimental PLF and TLF fragments, detected in the same ternary event are drawn as 
reference in the same plot. Data are compared with results of SMF+GEMINI calculation for 
Z=4 neck emitted fragments (hashed red histogram). Calculated data are corrected for the 
detector geometrical acceptance, thresholds and time resolution. The shape and position of the 
experimental data are well reproduced for all IMFs with charges from Z=3 to Z=8 (only charge 
Z=4 is shown in the plot). More details on calculations and comparison with experimental data 
can be found in [TT] . 




Figure 5. a) For the 124 Sn+ 64 Ni reaction, experimental < N/Z > distribution of IMFs 
as a function of charge Z for dynamically emitted particles (solid circles) and statistically 
emitted particles (solid squares); b) solid circles: same experimental data of Fig. 5a) for 
dynamically emitted particles. Blue hatched area: SMF-GEMINI calculations for dynamically 
emitted particles and asy-stiff parametrization; magenta hatched area: asy-soft parametrization 
for dynamically emitted particles, c) solid circles: experimental < N/Z > as a function of 
cos(6pnox) f° r charges 5 < Z < 8; empty circles: SMF calculation for primary fragment (asy- 
stiff parametrization); SMF-GEMINI calculations are indicated by blue-hatched area (asy-stiff 
parametrization) and magenta hatched area (asy-soft) respectively. 
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Figure 6. For the 124 Sn+ 64 Ni reaction, black empty histogram: velocity spectra Vji in the cm. 
system of dynamically emitted fragments with charge Z=4. Hashed red histogram: calculated 
velocity spectra (SMF+GEMINI) for the same reaction. Calculated data are filtered for detector 
acceptance. Experimental velocity distributions (yellow histograms) of PLF and TLF residues, 
detected in the same ternary event are plotted as reference. 

5. Summary and outlook 

We have studied with the 4tt detector CHIMERA the two reactions 64 ' 58 Ni+ 124 ' 112 Sn and 
I24,ii2g n _|_ 64,58]^j ^ ^ e same energy of relative motion (35 A MeV ). We defined a method to 
disentangle sequentially from dynamically emitted particles at midrapidity and to correlate the 
isotopic composition of intermediate mass fragments with their emission timescale. Dynamically 
emitted IMF shows larger values of <N/Z> isospin asymmetry and stronger angular anisotropics 
supporting the concept of "isospin migration" in neck fragmentation mechanism. Comparing 
the data for the reverse kinematics reaction with a stochastic mean field simulation, valuable 
constraints on the symmetry energy term of nuclear EOS at sub-saturation densities is obtained. 
A stiff E sym (p) behavior, with L «80 MeV corresponding to a linear density dependence of the 
symmetry energy, better reproduces our data. All these aspects open new perspectives for 
reaction studies with exotic beams. Our first outlook is in fact the possibility to plan new 
experiments using a 30 A MeV 68 iVi beam recently produced at LNS |23j . As a second outlook, 
we have recently proposed to study the 124 Xe+ 64 Zn reaction as compared with 124 Sn+ 64 Ni 
system where only the N/Z changes for the two systems with the same masses. We hope this 
study will permit to disentangle mass from isospin asymmetry effects evidencing the effective 
role of symmetry energy in the dynamics of the reactions. 
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